Robust Level Coincidences in the Subband Structure of Quasi 2D Systems 
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Recently, level crossings in the energy bands of crystals have been identified as a key signature for 
topological phase transitions. Using realistic models we show that the parameter space controlling 
the occurrence of level coincidences in energy bands has a much richer structure than anticipated 
previously. In particular, we identify robust level coincidences that cannot be removed by a small 
perturbation of the Hamiltonian compatible with the crystal symmetry. Different topological phases 
that are insulating in the bulk are then separated by a gapless (metallic) phase. We consider 
HgTe/CdTe quantum wells as a specific example. 
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Recently level crossings in the energy bands of crys- 
tals have become a subject of significant interest as they 
represent a key signature for topological phase transi- 
tions induced, e.g., by tuning the composition of an alloy 
or the thickness of a quasi two-dimensional (2D) system 
[2-|i|- For example, it wasproposed [f| and soon after 
confirmed experimentally [y, LD] that HgTe/CdTe quan- 
tum wells (QWs) show a phase transition from spin Hall 
insulator to a quantum spin Hall regime when the lowest 
electron-like and the highest hole-like subbands cross at 
a critical QW width of ~ 65 A; see also @, 0- 11 1. Here 
we present a systematic study of level crossings and anti- 
crossings in the subband structure of quasi 2D systems. 
We show that the parameter space characterizing level 
crossings has a much richer structure than previously an- 
ticipated. In particular, we present examples for robust 
level coincidences that are preserved while the system pa- 
rameters are varied within a finite range. Similar to the 
topological phase transitions characterizing the quantum 
Hall effect [Hj], the insulating Zi topological phases [l| 
thus get separated by a gapless (metallic) phase. Such 
an additional phase was previously predicted in Ref . [HI . 
Yet it was found that this phase could occur only in 3D, 
but not in 2D. Also, it was not clear which systems would 
realize such a phase. Here we take HgTe/CdTe QWs as a 
realistic example, though many results are relevant also 
for other quasi 2D systems 

Level crossings were studied already in the early days 
of quantum mechanics 14- 1(|. They occur, e.g., when 



atoms are placed in magnetic fields in the transition 
region between the weak-field Zeeman effect and the 
high-field Paschen-Back effect. Also, they occur when 
molecules and solids are formed from isolated atoms. 
Hund [13] pointed out that adiabatic changes of ID 
systems — unlike multi-dimensional systems — cannot give 
rise to level crossings. Von Neumann and Wigner 15 1 



quantified how many parameters need to be varied for a 
level crossing. While levels of different symmetries (i.e., 
levels transforming according to different irreducible rep- 
resentations, IRs) may cross when a single parameter is 



varied, to achieve a level crossing among two levels of the 
same symmetry, it is in general necessary to vary three 
(two) independent parameters if the underlying eigen- 
value problem is Hermitian (orthogonal). Subsequently, 
this problem was revisited by Herring [16] who found that 
the analysis by von Neumann and Wigner was not easily 
transferable to energy bands in a crystal due to the sym- 
metry of the crystal potential. Similar to energy levels 
in finite systems, levels may coincide in periodic crystals 
if the levels have different symmetries. Of course, unless 
the crystal is invariant under inversion, this can occur 
only for high-symmetry lines or planes in the Brillouin 
zone (BZ), where the group of the wave vector is different 
from the trivial group C\ . If at one end point ki of a line 
of symmetry a band with symmetry Ti is higher in energy 
than the band with symmetry Tj , while at the other end 
point k 2 the order of Ti and Tj is reversed, these levels 
cross somewhere in between ki and k 2 . Herring classified 
a level crossing as "vanishingly improbable" if it disap- 
peared upon an infinitesimal perturbation of the crystal 
potential compatible with all crystal symmetries. In that 
sense, a level coincidence at a high-symmetry point of 
the BZ such as the T point k — becomes vanishingly 
improbable. For energy levels with the same symmetry, 
Herring derived several theorems characterizing the con- 
ditions under which level crossings may occur. In partic- 
ular, he found that in the absence of inversion symme- 
try level crossings that are not vanishingly improbable 
may occur for isolated points k such that these crossings 
cannot be destroyed by an infinitesimal change in the 
crystal potential, but they occur at some point near k. 
Here we identify several examples for such robust level 
coincidences. This illustrates that level coincidences in 
energy bands can be qualitatively different from level co- 
incidences in other systems jl5| . 

Recently, several studies focusing on topological phase 
transitions recognized the importance of symmetry for 
level crossings in energy bands 0, |jl-Hcj|. Murakami et 
al. studied the phase transition separating spin Hall 
insulators from the quantum spin Hall regime, focusing 
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on generic low-symmetry configurations with and with- 
out inversion symmetry. They found that without inver- 
sion symmetry the phase transition is accompanied by 
a gap closing at points k that are not high-symmetry 
points. In inversion symmetric systems the gap closes 
only at points k = G/2 where G is a reciprocal lat- 
tice vector. Here we show that level crossings in quasi 
2D systems can be characterized by a multitude of sce- 
narios, taking HgTe/CdTe quantum wells as a specific 
example for which it is known that the lowest electron- 
like and the highest hole-like subbands (anti)cross for a 
critical QW width of about 65 A 0-0, . In most semi- 
conductors with a zinc blende structure (point group Td) 
the s-antibonding orbitals form the conduction band (IR 
Tg of Td), whereas the p-bonding orbitals form the va- 
lence band (Tg and T7 of Td). The curvature of the Tq 
band is thus positive whereas it is negative for the Tg, 
band. For finite k, the fourfold degenerate Tg states (ef- 
fective spin j = 3/2) split into so-called heavy hole (HH, 
m z = ±3/2) and light hole (LH, m z — ±1/2) branches. 
In HgTe, the order of the T 8 and T 6 bands is reversed: T 6 
is located below Tg and it has a negative (hole-like) cur- 
vature, whereas Tg splits into an electron (m z — ±1/2) 
and a hole (m z = ±3/2) branch [Hj]. HgTe and CdTe 
can be combined to form a ternary alloy Hg^Cdi-^Te, 
where the fundamental gap E$ between the Y§ and Lg 
bands can be tuned continuously from Eq = +1.6 eV in 
CdTe to Eq = —0.3 eV in HgTe with a gapless mate- 
rial for 1 w 0.16 [l8j ]. Tuning the material composition 
x thus allows one to overcome Herring's conclusion [l6| 
that a degeneracy at k = between two levels of different 
symmetries is, in general, vanishingly improbable. 

Layers of HgTe and CdTe can also be grown epitaxi- 
ally on top of each other to form QWs. At the interface 
the corresponding states need to be matched appropri- 
ately. The opposite signs of the effective mass inside and 
outside the well result in eigenstates localized at the in- 
terfaces We calculate these eigenstates as well as 
the corresponding subband dispersion E a (k) using a re- 
alistic 8x8 multiband Hamiltonian TL for the bulk bands 
Tq, Lg, and L7, which fully takes into account important 
details of E a (k) such as anisotropy, nonparabolicity, HH- 
LH coupling, and spin-orbit coupling both due to bulk 
inversion asymmetry (BIA) of the zinc blende structure 
of HgTe and CdTe as well as structure inversion asym- 
metry (SIA) of the confining potential V(z). For details 
concerning TL and its numerical solution see Refs. [U ED]. 
In the following k = (k x , k y ) denotes the 2D wave vector. 

The symmetry group Q of a QW and thus the allowed 
level crossings depend on the crystallographic orientation 
of the surface used to grow a QW [a (001) surface being 
the most common in experiments]. It also depends on 
whether we have a system without or with BIA and/or 
SIA. The resulting point groups are summarized in Ta- 
ble |U We show below that these different groups give 
rise to a rich parameter space for the occurrence of level 



TABLE I: The point group of a QW for different growth di- 
rections starting from a bulk semiconductor with diamond 
structure (point group Oh) or zinc blende structure (point 
group Td) for a system without ("sym.") or with ("asym.") 
SIA. 
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FIG. 1: (Color online) Subband states in a symmetric 
HgTe/CdTe quantum well (for k = 0) as a function of well 
width w calculated with an 8 x 8 Hamiltonian (a) neglecting 
BIA (point group D 4 h) and (b) with BIA (D 2 d). States trans- 
forming according to L^ of D ih (^6 of D 2d ) are shown in red; 
states shown in black transform according to L^ of D 4 h (L7 
of Aw). 



coincidences. For a proper symmetry classification we 
project the eigenstates of TL onto the IRs of the respec- 
tive point group [22|. In the following, all IRs are labeled 
according to Koster et al. (23|. As s pin- orbit coupling 
plays a crucial role for BIA and SIA [20] as well as for 
topological phase transitions [H-ll], all IRs referred to in 
this work are double-group IRs. For comparison, Table U 
also lists the point groups if the prevalent axial (or spher- 
ical) approximation is used for TL. In this approximation, 
BIA is ignored and different surface orientations become 
indistinguishable. 

First we neglect the small terms in TL due to BIA so 
that the bulk Hamiltonian has the point group Oh- In 
the absence of SIA, a quasi 2D system grown on a (001) 
surface has the point group D^h (which includes inver- 
sion) and all electron and hole states throughout the BZ 
are two-fold degenerate |22|. Subband edges k = in a 
HgTe/CdTe QW as a function of well width w are shown 
in Fig. HJa) . The HH states transform according to Tf of 
D^h. The electron-like and LH-like subbands transform 
according to Ly. As expected, the T 6 and Ly subbands 
may cross as a function of w. 

In the presence of SIA we cannot classify the eigen- 
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states anymore according to their behavior under parity. 
Without BIA the point group becomes C 4l) . HH states 
transform according to T 6 of C 4l , and electron- and LH- 
like states transform according to Tf, The level cros sing s 
depicted in Fig. QJa) remain allowed in this case 0, 24 1 . 



The situation changes when taking into account BIA. 
Without SIA the point group becomes D 2 d- In this case, 
all subbands transform alternately according to the IRs 
Tq and T7 of Did-, irrespective of the dominant spinor 
components. In particular, both the highest HH state 
and the lowest conduction band state transform accord- 
ing to T 6 of D 2 d so that around w ~ 65 A we obtain an 
anticrossing between these levels of about 2.9 meV (for 
k = 0), see Fig.QJb) With both BIA and SIA the 

point group becomes C 2v ■ Now we have only one double 
group IR IV Thus it follows readily that all subbands 
anticross as a function of a continuous parameter such as 
the well width. 

While BIA opens a gap at k = 0, level coincidences 
remain possible for some k ^ when the well width w is 
tuned to a critical value w Considering a (001) 
surface with BIA, we find, indeed, that for each direction 
(pofk— (k,4>), critical values w and k exist that give rise 
to a band crossing. Thus we get a line in k space where 
the bands cross when w is varied within some finite range. 
This result holds for QWs on a (001) surface with BIA, 
without and with SIA (as studied experimentally in Refs. 
0,0])- As an example, Fig.|2Ja) shows k in the presence 
of a perpendicular electric field E z = 100 kV/cm. 

In general, three independent parameters must be 
tuned for a level coincidence in a quantum mechanical 
systems [H| if the underlying eigenvalue problem is Her- 
mitian. While the multiband Hamiltonian H used here 
[2(]j is likewise Hermitian (not orthogonal), only two in- 
dependent parameters (w and k = |k|) are necessary to 
achieve the level degeneracy. We have here an exam- 
ple for the robustness of band coincidences under per- 
turbations that was predicted by Herring [l6| to occur 
in systems without a center of inversion (in multiples of 
four). It shows that level coincidences in energy bands 
can behave qualitatively different from level coincidences 
in other quantum mechanical systems [H| . We note that 
the band coincidences found here are not protected by 
symmetry in the sense that — unlike the other cases dis- 
cussed above — the group of k is the trivial group C\ con- 
taining only the identity. 

The situation is different for quasi 2D systems grown 
on a (111) surface. In the absence of BIA and SIA, 
the point group is D 3 d- HH states at k = trans- 
form according to the complex conjugate IRs © 
or IY © Tg , where © indicates that these IRs must be 
combined due to time reversal symmetry. All other sub- 
band edges transform according to Tf. In the presence 
of BIA and/or SIA the point group becomes Cz v . Then 
HH states transform according to the complex conjugate 
IRs r$ ®Tq. Electron- like and LH-like states transform 
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FIG. 2: Critical wave vectors k that give rise to a level coinci- 
dence in a HgTe/CdTe QW (a) on a (001) surface taking into 
account BIA (b) on a (110) surface neglecting BIA. In both 
cases a perpendicular field £ z = 100 kV/cm was assumed. In 
(a) the level coincidence requires a well width w = 66.1 A 
for k || [110] and w = 66.3 A for k || [110]. In (b) we have 
w = 60.9 A for k || [001] and w = 60.7 A for k || [110]. 



TABLE II: Irreducible representations of quasi 2D states 
(k = 0) on a (001) and (111) surface, starting from a bulk 
semiconductor with diamond (point group 0%) or zinc blende 
(point group T d ) structure for a system without ("sym.") or 
with ("asym.") structure inversion asymmetry. 
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according to T 4 . Thus it follows that on a (111) surface 
the HH states always cross the other states at k = as a 
function of w [similar to Fig. HJa)]. The IRs for different 
geometries starting out from a (001) or (111) surface are 
summarized in Table HU 

Finally we consider quasi 2D states on a (110) surface. 
In the absence of BIA and SIA, the point group becomes 
£>2/i- Here, all subbands transform alternately according 
to and with the topmost HH-like subband being 
and the lowest electron- like subband being T§ . A 
level crossing as a function of w is thus again allowed at 
k = 0. In the presence of either BIA or SIA the sym- 
metry is reduced to C 2v . While the point group in both 
cases is the same [2^, we obtain a remarkable difference 
between these cases. With SIA the level crossing oc- 
curs for a line in k space, similar to the (001) surface, 
see Fig. [2jb). With BIA we obtain a level crossing only 
for k || [110] with k « 0.0012 A" 1 and w « 62.5 A, 
thus giving an example for the level crossings occurring 
for isolated points k / as discussed by Murakami et 
al. [2j. These examples illustrate that the occurrence of 
level crossings at either isolated points or along contin- 
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uous lines in paramete r sp ace is not simply related with 
the system symmetry 25] . In the presence of both BIA 
and SI A (group C s ) we have the same situation as with 
BIA only, i.e., adding SIA changes the values of k and w, 
but we keep k || [TlO]. 

In conclusion, we have shown that a rich parameter 
space characterizes the occurrence of level coincidences 
in the subband structure of quasi 2D systems. In partic- 
ular, we have identified level coincidences for wave vectors 
k =/= that cannot be removed by a small perturbation 
of the Hamiltonian compatible with the QW symmetry 
16] . Taking into account the full crystal symmetry of real 
materials is an important difference between the current 
analysis and previous work that considered only lattice 
periodicity, inversion and time reversal symmetry. The 
full set of symmetries imposes additional constraints on 
the band Hamiltonian beyond the torus topology of the 
BZ that reflects the translational symmetry. These ad- 
ditional constraints generally reduce the number of pa- 
rameters that are required to obtain level crossings [16| 
so that robust level coincidences can be achieved even in 
quasi 2D systems. As quasi 2D systems can be designed 
and manipulated in various ways not available in 3D this 
opens new avenues for both experimental and theoretical 
research of topologically nontrivial materials. 

As a specific example, we have considered HgTe/CdTe 
QWs, where a particular level crossing reflects a topolog- 
ical phase transition from spin Hall insulator to a quan- 
tum spin Hall regime 043- The robustness of the level 
coincidences found here implies that these phases, which 
are insulating in the bulk, are separated by a gapless 
phase similar to the metallic phases that separate the in- 
sulating quantum Hall phases 12J. While in HgTe/CdTe 
QWs the range of critical well widths w giving rise to the 
metallic phase is rather small (about 0.1 monolayers), we 
expect that future research will be able to identify mate- 
rials showing larger parameter ranges that can be probed 
more easily in experiments. We note that our symmetry- 
based classification of level crossings is independent of 
specific numerical values of the band structure parame- 
ters entering the Hamiltonian %. Indeed, our findings 
are directly applicable also to other quasi 2D systems 
made of bulk semiconductors with a zinc blende or dia- 
mond structure such as hole subbands in GaAs/AlGaAs 
and SiGe quantum wells. In general, the k • p coupling 
between the LH1 (Tj of D ih ) and HH2 (T 6 ) subbands 
gives rise to an electron-like dis per sion of the LH1 sub- 
band for small wave vectors k 26]. If these subbands 
become (nearly) degenerate at k = 0, the coupling be- 
tween these subbands becomes the dominant effect. This 
situation is described by the same effective Hamiltonian 
that characterizes the subspace consisting of the lowest 
electron and highest HH subband in a HgTe/CdTe QW 
It can be exploited if biaxial strain is used to tune 
the separation between the LH1 and HH2 subbands 
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